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INTRODUCTION 

Today the subcritical assembly is a convenient tool of 
the nuclear scientist rik! engineer* It hrs several inherent 
advantages which rkc it particularly useful. It la much 
cheaper to assemble and operate then is r full scale reaotor* 
fho amounts of materials used pro much leas. It is Inherently 
safe, no extensive shielding need be provided, and work, o^n 
be carried on with none of the physical restrictions neces- 
sary near o reactor* 

A natural uranium, graphite moderated subcriticcl assea- 
bly hes been constructed recently at Iowa State College. This 
assembly will be useful in at le r st two capacities* 71rrt, 
it will provide extensive instructional opportunities, and 
tV.e chriice to compare classroom calculations with actual 
experimental results. Secondly, it will provide a research 
facility for numerous thesis investigations by graduate 
students * 

Before other studies could be undertaken, it was neces- 
sary to develop safe sad efficient procedures as well as to 
deter- iae the basic operating ehrracterintics of this assem- 
bly. It was with this urpoae In mind that the present thesis 
was undertaken* 
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RKVISV OF 'Inf LITEFIATURF 

The historical role played by the subcritical acse .bly 
in achieving the first self sustaining nuclerr chain reaction 
hre been given by Smyth (1). It Is summarized below. 

On December 2, 1942, mankind initiated bis first self 
sustaining nuolerr c'r nin reaction at the University of Chi- 
cago (1, p. 70). ibis event marked the culmination of rnonths 
of concentrated effort by many scientists and engineers. It 
c-uld not have occurred as it did without the subcriticsl 
assembly lending the way. 

By June 1940 it had been determined that the beat pos- 
sibility for a self sustaining ehein re ctlon lay in the 
thermal neutron fission of using a heterogeneous mix- 

ture of graphite and natural uranium (1, p. 09). It was 
known that in the fission of U more than 2 fast neutrons 
were produo ed, on the average, for every neutron consumed 
in fission* ho’ ever, most of these spare neutrons v-ore 
absorbed by other materials in slowing down, or escaped alto- 
gether, and thus were loot before they in turn were able to 
react with other uranium nuclei to produce fission. T v e prob- 
lem are to errands the fissionable material and moderator in 
such a manner that enough of the neutrons were conserved to 
sustain a chain reaction. 

Two main avenues of approach to this probler >ere pos- 
sible. An elaborate series of investigations into the nuclear 
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properties of various materials wsa initiated-, and theoretical 
calculations wert begun b*>sed on these data. however, since 
speed was t’ e controlling factor at that time, the major 
effort was directed into the empirical approach, ond here 
t. e subcrlticsl assembly became the primary me^ns of Investi- 
gating the problem. 

The very smell amounts of materials available, end the 
extreme coots involved, made full scale experiments Impos- 
sible. The method used then hnc to construct a pile of about 
one fourth to one third the size believed necessary for 
criticality. 

Fro* the results of experiments conducted with such s 
eubcrltical assembly it was possible to predict the character- 
istics of e similar full scale reactor. 

Since such a subcrltlcal assembly could not sustain a 
chain reaction an external source of neutrons vac placed near 
the bottom of the structure* T^e flux distribution in such 
an assembly then decreased exponentially with distance from 
the external source- This fact led to the subcritical assem- 
bly being called the "exponential pile". 

Doctors Fermi and Bxilnrd had suggested the use of 
graphite as a moderator for a chain reaction. Further, it 
was they who developed the lattice structure, in which lumps 
of uranium would be placed at regular intervals in a matrix 
of graphite (1, p- £3). This lattice structure bed definite 
advantages over a uniform mixture of uranium and moderator. 
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In July 1941 the first eutcritlcal fssembly wan set up 
fit Columbia university. It ms a graphite cube 8 ft. on a 
side, containing appro x i sur tely 7 tons of uranium oxide dis- 
tributed in lumps throughout the graphite. An s: ternal 
radium-beryllium neutron source v s placed ne^r the bottom 
of the nsBoitbly . Similar structures lucre set up in September 
and October of 1941. The multiplication footer, k, and t v e 
infinite multiplication factor, ka,, were computed in each 
ccse. For the third '’S-embly, k ^ rerohed 0.87 (1, p. 42). 

Tie problem remaining was to get k ^ to reach or surpass 
1*00 by using ; urer materials and possibly by using different 
lattice arrangements. 

within the following yerr a series of experiments was 
conducted using subcrltical assemblies. In each experiment 
the uranituB-graptoi te ratio was varied, oxides of Improved 
purities were used, lattice spacing was varied, and various 
sizes end shapes of fuel lumps were tried, T^e ninth pile 
of this series was constructed in July 194£. The graphite 
used In t- is case was sufficiently free of Impurities that 
it absorbed 20 per cent less neutrons then the best previously 
available commercial grahit®. Calculations baaed on this 
assembly showed that for a reactor of infinite size, k « 

1.007. An infinite multiplication factor greater t^nn 1.00 
had been reached (1, p. 63), and the subcritical assembly 
had been the means of rohievlng this goal- 

Since 1942 the cufccriticnl assembly hns been used 
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extensively in developing new reactor typ' s and in studying 
sesoolf ted problems. At the Hanford Atonic Products " , per* , tlon 
over 100 subcrlticfl experiments have been condu ted (?, p. 
309) . Ihese experiments have had as their objective the find- 
ing of e better water-cooled graphite lattice for plutonium 
production and providing improved knowledge of lattice theory. 

Argonne National Laboratory has been conducting a aeries 
of fast exponential experiments ( ft, p. 342). Thin series of 
experiments wre designed to obtain fundamental information 
about the physics of dilute fast reactors- In connection 
with the Cchool of Nuclear Science and Engineering, Argonne 
National Laboratory also operates a graphite moderated natural 
uranium subcrltlcal assembly rr, well as a heavy water moder- 
ated natural uranium assembly. 

In order to aid in the design of liquid metal cooled 
reactors, North American Aviation has conducted a series of 
Bubcritlcel experiments with a simulated sodium coolant (2, 
p- 309) . ihese experiments included both natural and en- 
riched urenium fuel elements in a graphite moderator. North 
American Aviation has also conducted a series of sufccritlcal 
experiments to investigate the b-sic physics of a heavy water 
moderated lattice using various fuel enrichments (2, p. 268). 

Beg inn lag in 1944 end 1945 a series of suberitieal 
experiments was conducted at the 0 nit TUdge National Labora- 
tory. The purpose of these experl sente was to determine 
values of buckling and other related characteristics for 
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various diameters of natural uranium fuel rods in combination 
■with a series of different volume ratios of light 'ater 
moderator to uranium (2, p. 185) • At the present ti .e the 
Oak. idge lohool of -o ctor Technology hrs a grabhlte moder- 
ated natural uranium cubcritical assembly for use primarily 
as an instructional tool (3). 

A unique application of the cubcritical assembly was 
mode during the construction of the Brookhnven ’"crctor. 

During a two-month period the partially completed reactor 
vse used for a aeries of cubcritical measurements . The re- 
sults of these measurements were used to determine critical 
buckling and other lattice constants for the full scale 
reactor (2, p. 305). Later a series of subcrltlcal experi- 
ments vi r »8 conducted at Drookhsvon to investigate the char- 
acteristics of light water moderated reactors with slightly 
enriched fuel (2, p. 184). 
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fh© flux distribution In r subcriticnl assembly does* not. 
Bstisly the wav® eqaetljn for o critlcsl rc"Ctor- no wever, 
if the aysc-bly la quite large th® thermal neutron flux dis- 
tribution at a distance from boundaries end from the source 
can be represented quite well by the equation 



lion. Although this equation Is valid only for a hoasofenecuo 
aye tea, it may be used for a heterogeneous systr.t without 
serious error- There will be localised Irregulfrities, but 
Equation 1 will give the overall neutron distribution- 

Mater! rl buckling is determined from raeryurerents of the 
ther-el neutron flux distribution throughout the sufccrl tierl 
aesenbly. iho critical else of r reactor of the sax® cokj^rJ- 
ti-n snd geometric configuration ten then be determined by 
equating &$®.*Jttrlcsl buckling to the isetoriel buckling deter- 
mined above- • 

Eor the usual rectangular oubcrltical assembly Equation 
1 is expressed as 




t 





rnd z are considered to te separable vpri^bler, 




then 



# r - 4*j •*- i 1 44- • «- « - . 

v a . 

*' ' •••-<- 
*•» a •*• ► m ••*![» » a • *,% 

+ H H 

[. ." ] •*► 



fend Equation £ tec ores 

1 

i dy ^ 



1 d 2 X + 1 d£y r 1 a£z 



X 



•r r; • o 



( 4 ) 



C X 



** (} y ^ 



Let tin,-. 



** *cy ri 

*■ M* - - & 2 
* * 

1 1 £{ « . a* 

p 



a d 2 : 



Y 



( 5 ) 

( 6 ) 

(?) 



s 2 



(S) 



-4 £ 2 , 

Zcu-tioa 4 cmj be written 

<* ♦ p - r * Bt. 

If the origin of the coor&in* te system, la vgaMSs-i to be 
at the center of t v, e Fide of the «<j»erabiy nr*ra«t t v e "xternrl 
neutron source, t^e bound? ry con it ion a een be expressed nc; 

1. the thcra&l neutron flux duo to V'Q ester;.:] scarce 
of fast neutrons etude • constant, jo> t;b# 

entire plena z ■■* 0, end 

!• Thernsl neutron flux (jP* 0 f t x a f p/2, f * | b/£, 
and £ * c. (a, b, ana c *»re overall dlnenai one of 
the aaseth.ly, including extra* olrtion diatrncs.) 
for any ^iven value of s, t* e flux lo n srxiauiB at y * 0, y n 
0, and is ayawt® trie oily dl tributed. It then follows th«t 
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( 12 ) 



Substituting in Equation & leads to t^e form 





(13) 



The solution for Z is of the fora Z( z ) * F slnhtf(c - z) , 
end the general solution nay be expressed as 



Combining Equations 9, 10 end 14, the general solution of 
the wave equation becomes 



where A D , C R , and F Cil are combined into A^ n * 

If the external neutron source Is ot least two slowing 
down lengths below the z » 0 plane, then it can be assumed 
that virtually all neutrons have been thermalized by the 
ti^e they enter the aubcritical assembly. Under those con- 
ditions the only term of Equation 15 thrt is of any signify 
cence is the fundamental mode (m » 1, n * 1). Equation 15 
then reduces to 



z cn “ ^ an ninh ^^(c - z) 



(14) 
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This can be further reduced to 
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<P . F" e * c je~ t - o" r ° e*^ (18) 

for given physical dimensions, e is n constant, and 
Equation 18 becomes 

<P =. c' e-T^ [l - e- a5f(0 -^ (IS) 

If the vertical dimension, c, is very large, ?nd if z is 
not allowed to approach c, then £]. - e~‘ * -zjj approaches 
unity. Equation 19 simplifies to 

^U) - c' o" rz (0 0 } 

Equation CO indicates that the thermal neutron flux distribu- 
tion decreases exponentially tfetth increasing z. Thus s plot 
of distance t versus the natural logarithm of the flux (or 
of a quantity proportional to flux) should be a straight 
line of elope V . Near the top of the assembly the points 
will deviate from a straight line due to neglecting the end 
correction terra jj. - e “2 ? ( c~z*)J ln qq Ue t;i on oq. 

The slope, Q , is really cn ^ Equation 13 can be 

written 

(JL) + (*L) 2 - tf 2 * >2 ( 21 ) 

a 

This equation can be solved for arterial buckling, all 
other quantities now being known. By netting this value of 
material buckling equal to geometrical buckling, the critical 
size of the full scale counterpart of t> e subcriticnl assem - 
bly can be determined. 
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DKUCRIFflOJI CP APPARATUS A ,D EXPERIMENTAL PROCEDURE 

Figure 1 Is an overall view of the eubcritlcnl asse bly . 
Figure £ shows Important dimensions and. the arrangement and 
numbering of the foil slots. The assembly consists of 
cylinders of reactor grcde graphite 60 In. long and 7 in. In 
diameter. These cylinders hfd t^e sides cut flat to form 
squares 6 in. ecross with rounded corners. There nre stacked 
in a square lattice, ten rows wide end fourteen rows high. 

Tie top five rows rre only & in. high by 6 in. wide, making 
t e overall dimensions of the assembly 60 by 60 by 79 In. 
high . 

This pile is set on a wooden base, and an external 
neutron source is centered under it. 'ater filled aluminum 
tanks entirely surround the external neutron source, provid- 
ing adequate shielding against both gamma and neutron radia- 
tion. 

The rounded corners between the graphite blocks present 
117 holes approximately 1 1/2 in. across for insertion of 
fuel elements. The spacing of these holes allows a Mittimus 
fuel lattice of 6 in. when all holes are filled. By remov- 
ing alternate fuel elements, or entire rows of fuel elements, 
the lattice spacing may bo varied from 6 in. to 8.5 in. 

(Q.484 in.), 12 in. or larger if desired. 

Thirty-one small horizontal slots piorce the assembly 
from the east face (Figure 2). T^ese slot^ provide for the 
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1. Subcritical assembly 
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X INCHES 



Figure 2* Subcritlcfll assembly - dimensions end 
foil holder arrangement 
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Irradiation of indium foils to use in determining the neutron 
flux at various points throughout tho pile* One vertical hole 
is available for possible simulation of control rod effects- 
All vertical sides and the top of the assembly ere covered 
with 0.010 in* sheet cadmium to present essentially a “black 
boundary 8 to neutrons. 

The indium foils used in this experiment are 1.0' in. by 
1.50 in. by 0.003 in. thick. They ore mounted on aluminum 
plauchets with DuPont “Duco 8 cement, each with the weight and 
foil number inscribed on the bock. The average weight of 
the 160 foils in use is 0*5943 grams. Aluminum foil holders 
allow placing of the foils in the previously described slots 
at 3 in. center to center spacing rcross t’^e assembly. For 
measuring feat neutron flux indium foils con be wrapped in 
0.010 in. cadmium sheet. 

Foil counting apparatus consists of a suitable decimal 
scalar and either a mice end window Geiger-Hueller counter or 
a thin glass wall Gelger-Kueller counter. 

The use of indium folia to measure neutron fluxes makes 
use of the so-called activation method. In this method the 
indium is exposed in the neutron flux. Indium hre a high 
cross section for thermal neutrons end becomes radioactive, 
emitting e beta particle with a 54 min. half life. The 
neutrons nre quite difficult to me. sure directly, but t^e 
beta activity of the indium la proportional to the neutron 
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activity, end beta particles can be detected and counted 
quite readily. 

For each, experimental run indium foils were placed in 
the assembly and Irradiated for a minimum of 5.4 hr. or six 
half lives; in most cases the lrradirtion time was much longer 
than this. This procedure insured an activity of at le'-st 
96.4 per cent of noxicium of the 54 min. activity. Foils vere 
then removed individually or in groups of not more than 
eight, depending upon the particular experiment. The foils 
were counted in a nearby room having a much lover background 
than the room containing the cubcriticel assembly. Each foil 
was counted for 3 r.,in. Then these counts, less background, 
v.’ere corrected Dock to the tire of removal from the assembly 
to give net activity in counts per minute- In each instance 
In which more than -one foil was removed for counting at a 
given time the foil expected to have the least activity was 
counted first. This helped to minimize distortion of low 
counting rates by large corrections to time of removal* 

No foil was counted sooner than 3 min. after removal from 
t' e assembly. This allowed s 13 sec- indium activity to 
decay to lose than 0*01 per cent, end only the 54 .in. beta 
activity remained. 

The majority of the data in this set of experiments was 
obtained using a 100 millicurle radium-beryllium external 
neutron source. This was the largest source available at 
the time the experiments were begun. Five, 1 curie 
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plutonlum-cery Ilium sources became available lrtrr, and two 
sets of crta were obtalnod using this larger neutron source. 

The first experimental arrangement consisted of V'e 6 in. 
lattice spocixig with the 100 mllllcurie neutron source- A 
complete survey of the lower horizontal row of foils • lun the 
entire vertical row of foils wrs texen in this configuration, 
foils were spaced at 6 in. intervals from the east face to 
t v e centerline of the assembly. Next s vertical survey vsm 
nade usin, cadmium covered center line foils to determine 
Ue cadmium ratio. The third set of data was obtained with 
the same G in. lattice, but with only the £2 center line 
foils loaded. This was to Investigate the extent of flux 
depression caused by the 13£ foils of t^e original loading- 

The lattice spacing was then changed to 8-5 in. rnd 
with oil foils loaded only the 18 in. and -30 in- (center lina) 
foils were counted. A set of foils was also placed conuletely 
through the assembly in the y direction. Date obtained from 
these folia were used In determining any flux variation due 
to the unsyxmr.otrlcal foil lording employed. 

Fuel loading was changed again, and the 18 in. and 30 
in. flux measurements were repeated with a 12 in. lattice 
spacing . 

The 100 mllllcurie radium- beryllium. neutron source was 
replaced by a cluster of five, 1 curie plutonium- beryllium 
neutron sources. f*o further experimental runs were con- 
ducted. First a survey w r e taken wit^ all fuel elements 
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renoveu. These data were used in determining the effect of 
subcrltical atoltl plication • Final ly, r set of a at a wes col- 
lect e< using the 8-5 in* lattice with the larger neutron 
source. These data provided s correlation for the 8. A in. 
lattice between characteristics obtained the aaall 

external neutron source and chsrsctaris ties obtained v/lth 
t^e large oxterr.nl neutron source. Again one row of folia 
was extended entirely through the assembly to investigate t v c 
#ya»etry in the y direction. 

l v- a indiuu foils were positioned in a holder ^trroxi- 
m* tely C.^.0 in. fror the end window Geiger- .roller tube in a 
constant position relative to t^e tube. In each c^ee the 
sr«a of t^e foil was greater t^*n V'& '’rep subtended by the 
end window tube, and no correction was «»de for weight v*ri- 
»tione due to the slight area variations among th© foils. 

A otrontlua-y ttrlir. aanple was used as n atari era source, 
end this a&aple rrs counted r, t least once during f;n c* run. 

This procedure wen followed to insure th p t no change in the 
experimental armuguuent or deviation, in the Geiger- cellar 
tube characteristics would creep in undetec ted end invalidate 
the correlation between data obtained on different occfslor.e. 
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DISCUSSION OF RESULTS 

The results of this series of flux measurenento are pre- 
sented gx’cphically so Figures 3 through 14. The data on 
which these figures ere based are tabulated in the Appendix. 

In a number of the figures a horizontal flux distribution 
is plotted. Since the theoretical flux distribution bro the 
sh pe of r cosine curve, cosine curves were faired through 
th<- experimental points rather than attempting to fit some 
other smooth curve to the points. 

A detailed survey was made with the 6 in. lattice and 
the small external neutron source. Due to the small activity 
of this source, counting rates were quite low and standard 
deviations were necessarily large. Extended counting period* 
would have lowered deviations, but this was not practical 
due to the excessive decay of the 54 sin. beta activity while 
counting over a longer period. Figure 3 displays the flux 
distribution in the x direction as compared to the theoretical 
cosine distribution. Vertical height of the symbols indicates 
standard deviation. It la noted that the scatter of experi- 
mental points decreases end the flux distribution conforms 
more closely to the plotted cosine curves as y increases 
toward the center line of the assembly. 

It was desired to determine whether o large number of 
indium foils present in the assembly produced a mersunble 
effect on the overall flux level. Figure 4 shows that any 
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figure 3. Horizontal flux curvey, 6 inch lattice, 
2 « 30 inches 
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Figure 4. Horizontal Flux survey for light rrv3 heevy 
iridium loods, y « 30 inches, z a 30 Inches 
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variation due to foil lording in within the 11 Its of normal 
scattering ., and only one cosine curve is plotted throug’ both 
sets of experimental points. Ibis conclusion is also sup- 
ported oy the drtjn of Figure &. Figure 6 Illustrates the 
syra .etry In t ' t y direction oven though '.any -ore fol3.s ’-ere 
pieced in the e*«t hall' of the asseUsly t' an in the vest half. 

In order to determine the degree to v v ich the fast source 
neutrons had been therriollsod, a ver-ticnl flux survey was 
ascie using cndralun covered iridium folia. Figure 0 displays 
these data contrasted vlth the coap retie <2.«Jio obtained with 
bare indium folia- fi constant, ratio of t'eco two fluxes 
•would indicate that all of tbs fast source neutrons hsg been 
thermal ized (4, r. l£4) . Ihc date indicate th-t t>st neutrons 
still re Eta in for none distance ur fro® the fcottow of t*>e 
assembly, i v -e apparently changing ratio nerr t*’« top r ? \y be 
due to statistical scattering of the points because of l h e 
very low counting rates involved. It ie believed that further 
investigation using the large external neutron source would 
show this trend to be in orror. 

Figure S comperes the flux obtained with each of three- 
lattice spacing©. Ihe 8.5 in. lattice gives the highest flux 
levels of the three lattices investigated. Home of th* Arts 
ottelnefi *Ith the 18 in. lattice is believed to be In error, 
as t^ese were the la t counts reads before an equipment failure. 
Since the nature of the fc-lluro wan such that the counting 
rote became extremely high Just before complete failure, 
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Figure 5. Cadrcluu ratio - 6 Inch lattice, x » 27 inches, 
y * 30 inches 
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these experimental points rre believed to indicate rn un- 
reasonably high flux, Therefore, the experimental points at 
x * 2? In. sad greater were disregarded. Instead, the cosine 
curve was fitted to points for x less than 27 in., and the 
curve drawn symmetrically . The dashed line indicates the 
extrapolated portion of the curve. 

frior to the completion of all experimental determ* na- 
tions a lrrger neutron souroe became available, ftor.e data 
were obtained using this larger source with t v e 8.G in. 
lattice. These results are presented both to confirm the 
findings obtained with the small neutron source end to demon- 
strate t v e great increase in flux levels obtained by use of 
a larger neutron source. 

Figures 7 and 8 display the flux distribution in the x 
criu y directions, respectively, with each of the neutron 
sources in piece. In both cases the higher flux levels, 
corresponding to the larger neutron source, show little 
scatter of experimental points, excellent symmetry, and o 
good conformance to the theoretical cosine distribution* As 
has been noted previously. Figure 8 Indicates no measurable 
deviation from the cosine flux distribution due to unsym- 
metrlcal indium foil loading in the y direction. 

One set of data was obtained with all fuel elements 
removed from the assembly. The flux levels present under 
these conditions, when coicprred with the flux levels present 
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Figure ?* Horiaontel flux survey, x direction, 6.6 inch 
lattice, different neutron sources 
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Figure 6. horizontal flux survey, y direction, 8.5 Inch 
lattice, different neutron courcen 
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in the 8.5 in. lattice, indicate the effect of subcriticnl 
multiplication* Figure 9 shows that within 30 in. vertically 
of tie base of the assembly the empty graphite lattice h^s n 
higher flux level than does the one with fuel elements in 
place* This is an indication that more neutrons are being 
absorbed by tht fuel than rr.e being produced by fission within 
the fuel. Above -30 In., however, the number of neutrons 
produced by fission in the uranium exceeds the number absorbed 
in the fuel, and the flux level with fuel elements in plsae is 
higher than the flux level of the empty assembly. Figure 10 
compares the flux level at the 30 In. plane for various 
values of x horizontally. Ae stated above, the overall 
effect of the fuel elements is essentially zero at this dis- 
tance from tho base. 

Ons of the uoot important characteristics which can be 
determined experimentally by mesas of a subcriticnl assembly 
is critical buckling. As developed previously, when flux (or 
ft vrlu. proportional to flux) is plotted on a logarithmic 
scale versus distance from the external neutron source on a 
linear scale the slope of the line thus obtained is gamma. 

If the physical dimensions of the subcrltical assembly ere 
known, critical buckling can then be determined by using 
Equation iil . 

In theory, critical buckling can be determined by first 
computing the infinite multiplier tion factor, k^, and then 
solving the transcendental equation 
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Figure 9* Effect of euboriticnl multiplication, B.5 
lattice, x «* £7 Inches, y « -50 inohoo 
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Figure 10. Effect of subcrltionl multiplication, 6.5 inch 
lattice, y « <30 Inches?, a « -30 inched 
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(22) 



The Infinite multiplier tion fector is computed from the 
four fector formula 



gives a value for natural uranium fuel of 'h * 1-553. Cross 
sections used are those given by Kaplan (5, j» . 490). From 
e figure given by Glass tone (6, p. 195), the fast fission 
factor, 6 , for the size fuel elements used in the subcritlcai 
assembly is 1-030. These two factor® are independent of 
lattice spacing and remain the s f ne for all three lattices- 
The resonance escape probability, p, varies with the 
relative amounts of fuel and moderator present, and It is 
therefore different for e ' ch lattice spacing. L^e method of 
Glnsstone (6, pp- 191-192) was used to determine p- 'eights, 
volumes, and densities of fuel, cladding and graphite were 
determined from representative samples. The density of the 
graphite was computed to be 1-56 g/ca 3 - The overall dimen- 
sions of © fuel element were measured to be 8-4u in. long by 
1-080 in. in diameter- Each fuel element weighed £-010 kg, 
giving it an overage density of 18-04 g/c* 3 • Each fuel 
element coiitained a rod of natural uranium 8.00 in. lonf rnd 
I- 00 in. in diameter* This left cladding of 0-040 in. thick 
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Co-puting 'Tj from the relation 
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with b • £00 in. end pi? tea. r-t density of the "luninuro 
cladding end binder (assumed the '•ame) ' rs taken to be that 
of alu. inuu, 2.71 z/ T v e density of the uranium was then 

dC 

computed to be 19-0 g/ca - ' • 

Using these data, t v e nqu«re heterogeneous lattice v-*s 
treated rs a system of equivalent c lindrl c*>l unit cells (4, 
p. , and p w^g detrained for e*>oh case. 

fbe fourth factor, fieri?,' utilisation, f, •’Iso varies 
wit! lattice rpecinfc. Thermal utilization w*s oa-euted *»s 
si own by iurrey (7, p. 37} , end & l «in the -qa^re lattice was 
treated * 6 « system of equivalent cylindrical unit cello. 
Allowance was made for th e presence of fuel, cladding, «ir 
spree, ciid graphite aoderctor . Tor l' e size fuel element 
used, the disadvantage factor, f, was covuted to be 1*0991. 

CoEputod values of ->| , 6 , ,, f, •Rd * ore tabulated 
in lable 1- 



iable 1. ibiorstical calculation of critical buckling 



Value 


6 in. lattice 


6.5 in. lattice 


12 in. lattice 




1.3W 


1.353 


1 . 353 


6 


1.030 


1.030 


1.030 


p 


0.808 


0.900 


0.94 7 


f 


0 .928 


0.873 


0.737 


k oo 


1-047 


1.098 


1.015 


B^( cc“’-) 


83 x 10" 6 


13? x 10“ 6 


15.5 x 10“ 
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having determined k^ , n first approximation for critical 
buckling con be obtained from Formula 25 w v< ieh Is v-lid for 
a large reactor 

8 2 « a K.r .. l (25) (4, p. 216) 

l 2 

r-w 

Feral age, (_ , for graphite Is 550 on • For a heterogeneous 
lattice, age may be taken os 10 per cent greater than that 
for the pure moderator ( 6 , p. 2C1), giving a value of *385 cm c . 

O 

Using values of f Just computed, L*“ was obtained from the 
formula 

L 2 * L 2 (l - f) (26) (4, p. 260) 

Values of buckling obtained by Equation 25 were then sub- 
stituted in the transcendental equation, Equation 22* The 
maximum error wes 0.3 per cent- Blight adjustments gave 
values of buckling correct within 0.1 per cent. 

The results of the above calculations for all three 
lattices ore tabulated as Table 1. 

The results of experiments! determination are presented 
as Figures 11, 12, end 13. The buekllngs computed from these 
figures ore listed in Table 2. 

In the interest of clarity, ell buckling data Vve been 
displayed on a common figure. Figure 14 compares theoretical 
buckling from Table 1 with experimental values deter ined 
from Figures 11 through 13, as tabulated in Table 2. 

Figure 11 compares vertical flux distribution for the 
three lattices at x « 27 In. snd y a 18 in. Ho attempt was 
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Figure XI. Lattice comparison - vertical, x * f7 inches, 
y * 18 inches 
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Figure 12. Lattice cocprrlson - vertical, x =» L7 Inchon, 
y ** 30 inches 
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Z INCHES 

figure 13. Vertical eurvey, 8.6 Inch lattice, different 
neutron sources 
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• Theory 
o Y = 1 8 inches 
A Y =30 inches 




t'l L ur« 14. Buckling veroua lattice spacing 
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Table 2* Experimental values of critical buckling 

“cm" 2 •> 



Configuration 




0 in. 
lattice 


8.5 in. 
lattice 


12 in. 
lattice 


Small neutron 
x » 27 in., y 


source, 

» 18 in- 


141 x 10” 6 


181 x 


io“ r> 


Hot 

determined 


Small neutron 
x » 2? in., y 


source, 

* 30 in. 


161 x 10“ 6 


65 x 


10-6 


19 x 10“ 6 


Large neutron 
x = 2? in., y 


source, 

* 30 in. 


Mot 

determined 


65 x 


10“ 6 


Hot 

determined 



made to estimate buckling l'or the 12 in. lattice r s the 
experimental points arc too badly scattered. This is be- 
lieved due to the equipment failure mentioned previously. 

Figure 12 presents the same comparisons cl y =» 30 in. 
This time the data for the 12 in. lattice appeared to be .more 
reliable, and buckling was computed for ell three lattices. 
Figure 13 presents the data for the 8.5 in. lattice from 
Figure 12, and compares it with data obtained for the ra7se 
lattice with the large external neutron source. Parallel 
straight lines have been plotted for the two nets of d«ta, 
and the experimental points fit these lines very well. It is 
noted that the slope of these lines gives © lower value of 
buckling for the 8.5 in. lattice at y = 30 in. than did the 
data of Figure 11 taken at y «* 18 in. 

It should be noted that a very snail change in the slope 
of a line makes a marked difference in the value of critical 
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buc&llng computed from that elope- For tMs reason, repeated 
experiments at a high flux level would be necessary before 
completely reliable values of critical buckling could be 
deter lined for a given lattice arrangement. The experimental 
critical buckllngs shown in Figure 14 ohould be considered to 
be indicative of general magnitude rather than highly accurate 
values. 
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C: i CLU3I0W0 

The operating chcrnoterletico of the eubcritlcrl assembly 
under Investigation agree quite ,- ell with the accepted theory, 
ire 6*6 in. lattice le the optimum lattice, resulting in 
higher flux levels end the largest value of eritlccl buckling. 
All experimental buckling determinations proved to be posi- 
tive, which indicates thet a full scale rerotor could be 
constructed of the came configuration rs either the 6 in., 
fl.5 in., or 12 in. lattice. There is no apparent measurable 
flux depression due to Indium foil loading, even when foils 
are loaded unsynsmetric^lly • Horizontal flux distributions 
ere symmetrical in both the x end y directions, end within 
the limitations imposed by statistics they agree well wit^ 
the cosine distributions* predicted by theory. 

The euberitiesi assembly could be operated satisfactorily 
using the smaller (100 nsillicurle) external neutron source. 
However, the scattering of the experimental points obtained 
with this source is quite severe, and statistical deviations 
are very large. With the larger neutron source in place, flux 
levels ore higher by a factor of five or ore, with corre- 
spondingly smaller deviations and less scattering of experi- 
mental points. 

The fast neutrons from the external source are not 
entirely theraallzed until they have penetrated approximately 
30 in. into t v e assembly. 
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This subcritical assembly, with the large external neutron 
source In place, can be used to great advantage os on instruc- 
tional tool* Although t v e experimental results obtained were 
not extremely accurate, it is believed thrt the dr ta can be 
reproduced with sufficient consistency and eccurrcy to cor- 
relate experimental results with classroom calculations* 
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SUGGESTIONS FOR FURTHER STUDY 

Further study utilising the sene subcrltlesl assembly 
used in this work, might well be carried out in one of three 
general directions- 

In the first of these categories lies the determination 
of various other properties of the assembly, such as diffu- 
sion length, Fermi age, and effective cross section- Methods 
of determining these quantities experimentally are given by 
Hurray (?, p. 105) - Hone of these properties h*r. been deter- 
mined experimentally for this particular aubcritlcal assembly. 

The second group of possible experiments includes all 
those in which some basic physloal parameter in varied and 
the effect on nuclear properties io observed. Experiments 
such as the study of local flux depr 33ione under various 
partial fuel lo^ds, simulation of coolants, temperature 
effects on nuelc°r properties, variation of neutron source 
geometry, end simulation of control rods all cone under t^ie 
classification. For any lattice spacing other than 6 in. 
there is an appreciable volume of air contained In the empty 
fuel channels. Flux cc fiureaents should be made with these 
fuel channels filled with graphite to determine th© effect, 
if tny, of the empty channels. 

The last type of experiment for v v ieh t v, e oubcrltlcal 
assembly night be used is one In which the assembly rets as 
a testing facility, rather than being investigated itself. 
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Thio might inolude feasibility studies of e new type of 
neutron counter, or the development of s movable, continuously 
monitoring ana recording survey Instrument. This also could 
include the experimental justification of newly developed 
theory. An example of this type of application is the devel- 
opment of some now pjrid less critical method of experimentally 
determining critical buckling. 
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Table 3. 


Key to 


tables In 


Appendix 








Neutron 


Lattice 




Table 


Run 


source 


©pacing 


Comments 


4 


1 


small 


6 in. 




5 


2 


small 


0 in. 


Cadmium covered 


6 


3 


f>»Rll 


6 in. 


Light Toil load 


7 


4 


snail 


i . v in . 




e 


5 


email 


12 in. 




9 


6 


large 


none 


All fuel re oved 


10 


7 


large 


8-5 in. 
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Table 4. Run 1 



Position 


Saturated 
cctivi ty(c/~a) 


Position 


Saturated 
sctivl ty( c/aj) 


0 


4 


P-5-30 


115 


A — &— 1> 


5 


G-5-0 


14 


A-5-1*. 


15 


G-6-6 


•37 


A- 5-18 


21 


(5-5-12 


64 


A- 5- -4 


25 


0-5-16 


73 


A- 5- 50 


26 


G-5-24 


91 


B-5-0 


0 


C— 5-30 


102 


t-6-0 


11 


H-5-0 


£ 


n-b-li. 


21 


H-5-6 


29 


8-5-16 


58 


n- 6-12 


43 


n- 5- 24 


49 


H-5-18 


57 


S- 5-30 


4 5 


j ;-5-24 


62 


C-5-0 


10 


H-5-30 


80 


C-5-6 


2-3 


1-5-0 


4 


C-5-12 


40 


1-5-8 


19 


G-5-18 


66 


1-5-12 


41 


C-5-^4 


75 


I -5-18 


42 


C-5-30 


89 


1-5-24 


47 


D- 5-0 


6 


I- 5-30 


54 


5-5-6 


25 


J-5-0 


l 


L-6-1L 


50 


J- 5-6 


6 


D-5-18 


73 


J-5-12 


14 


D-5-24 


89 


J-5-18 


16 


i>— o— 3w 


99 


J- 5-24 


25 


E- 5-0 


7 


J-5-30 


£3 


t-C-6 


38 


E-13-0 


V 


£- 5-11 


72 


F-13-6 


3 


E-6-18 


76 


E-13-12 


0 


1-5-24 


105 


E-13-18 


1 


E-5-30 


129 


t -13-24 


4 


F-5-u 


10 


-’-13-30 


4 


P-5-6 


24 


-12-0 


4 


P-5-12 


68 


-12-5 


6 


r-e-18 


80 


I -12-12 


6 


F-5-24 


91 


F-12-18 


5 
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Table 4. (Continued) 



Position 


Saturated 
activi ty( c/r») 


Position 


Fnturs ted 
activi ty( c/") 


F-12-C4 


6 


F.-6-18 


50 


E-U-50 


S 


7-6-24 


74 


£-11-0 


0 


£-6- 30 


86 


£-11-6 


0 


E-4-0 


14 


£-11-12 


2 


’ -4-6 


52 


E-.U-13 


10 


K-4-12 


87 


£-11- £.4 


8 


f-4-ie 


106 


i- 11-30 


8 


» -4-24 


170 


K-10-0 


1 


-4-50 


188 


£-10-6 


7 


P-3-0 


11 


£-lu-l£ 


16 


-3-6 


70 


-10-18 


17 


FI- 3- 12 


117 


£-10- 14 


19 


-3-18 


195 


t-10-30 


19 


r -3-£4 


256 


£-9-0 


0 


P-3-30 


269 


-3-6 


9 


£-2-0 


18 




*,i 


1-2-6 


71 


F.-9-18 


15 


2-2-12 


141 


E-9-C4 


20 


P-2-18 


235 


£-3-30 


25 


F- 2-24 


341 


£-8-0 


7 


£-2-30 


368 


E-3-ft 


10 


P-1-0 


18 


E-8-12 


c2 


-1-6 


88 


£-8-18 


* 2 


P-1-12 


167 


K-8-^4 


29 


£-1-18 


315 


'•-8-30 


58 


F -1-24 


446 


E-7-C 


3 


r -l— 30 


504 


E-7-G 


10 






E-7-1L 


£4 






£-7-18 


25 






E-7-24 


48 






£-7-30 


58 






£-6-0 


4 






E-0-6 


20 






£-0-10 


63 
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i able 5. Run 2 



i ocl tion 


Sr turated 
aetlvity( c/a) 


Fonltion 


2 o turs ted 
ectivity( c/ra) 


K— _.-30 


10 


F-3-30 


09 


R.-9-3Q 


12 


E-2-30 


150 


E- 10-30 


1 


2-1-30 


201 


r-11-30 


7 






E-l-SO 


0 






F.-l 5- 30 


4 






'.-7-30 


17 






2-0-30 


£4 






E-5-30 


31 






k-4-30 


51 







TsL-le 6. Run 3 




Sp turated 




Saturated 


Position 


cctivlty(c/o) 


Position 


sctlvi ty( c/e) 


£-1 >-50 


4 


6-5-30 


106 


E- 12-30 


7 


F-5-30 


103 


C-ll-30 


14 


6-5-30 


122 


- 10-30 


15 


-4-30 


176 


£-9-30 


24 


1 

* 

o 


2CQ 


£-6-30 


34 


E-2-30 


339 


E-7-50 


47 


’ -1-30 


525 


/ i -5-3 f J 


19 






R— 5— 


47 






C-5— 30 


82 






J-5-30 


16 






3-5-30 


46 






r.-5— 30 


70 






E-C-3C 


75 






D-5-30 


99 
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Tpble 7. Run 4 



Posi tion 


SRturrted 

actlvIty(o/r3) 


position 


Saturr ted 
activity (c/c) 


E,-l3-lb 


1 


G-5-18 


9.3 


-1-5-30 


3 


G-5-30 


10 3 


F-1S-18 




F-5-18 


108 


t -12-30 


£ 


F-6-30 


132 


e-ii-is 


6 


r-5-18 


100 


L-ll-5u 


15 


E-5-30 


139 


C-10-13 


1£ 


E-4-18 


152 


£-10-30 


15 


F-4-30 


£04 


L-9-18 


51 


C— 3-18 


£18 


1-9-30 


* >o 


r-3-30 


352 


E-e-16 


39 


r-2-is 


310 


-8-3C 


33 


0-5-30 


475 


E-7-16 


46 


6-1-18 


398 


\ -7-30 


66 


8-1-30 


610 


#.-8-18 


67 


0-5-0 


12 


E-e-30 


83 


7-5-6 


47 


/— 5-18 


34 


7-5-12 


74 


A- 5-30 


19 


F-5-24 


126 


8-5-18 


46 


7-5-36 


12S 


8-5-30 


69 


F-6-4S 


120 


C-S-18 


71 


7-5-4 9 


80 


C-5-30 


115 


7-5- 54 


45 


8-5-18 


05 


r-5-60 


20 


0-5—30 


159 






J -5-18 


.33 






J-5-30 


18 






1-5-18 


56 






1-5-30 


68 






..-5-18 


78 






/*— u— 30 


97 







52 



TrLlo 8. Hu:> 5 



Pool tlon 


Spturrted 
ac tivity( c/n) 


Position 


Sr turf tel 
activity ( c/m) 


r.-l 3-18 


8 


J-8-30 


IQ 


E-io- 38 


6 


1-5-18 


51 


-lx-18 


6 


1-5-30 


53 


l.-l^-30 


6 


8-5-10 


7C 


I.-11-18 


14 


0-5- 30 


102 


-lx- 30 


7 


G-5-18 


105 


-10-18 


10 


G-5-30 


185 


-18-30 


17 


F-5-18 


104 


•0-9-18 


18 


F-5-30 


143 


£-9-51 


£8 


F -5-18 


ICO 


-8-18 


28 


* -5-30 


135 


*-8-3u 


46 


E-4-18 


CO 3 


1-7-18 


33 


7.-4- 30 


?41 


E-7-30 


40 


~- 5-18 


geo 


B.-G-18 


62 


F-3-50 


374 


E — 6— 58 


89 


E-2-18 


340 


A-0-1G 


18 


E-C-30 


550 


A- 0-30 


19 


E-l-18 


4C6 


8-5-18 


48 


r-i-30 


663 


B-o-30 


53 






C-5-18 


70 






C-0-§0 


90 






5— 0— lb 


9c 






I. -0-30 


113 






J-5-18 


14 
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Table 9. Run 6 



Poel tion 


Saturated 
activity (c/e) 


Po nit ion 


Saturated 
nctlvity( c/m) 


E-13-3^ 


9 


J-5-30 


135 


E-Il-30 


22 


r-5-o 


70 


E- 11-50 


33 


F-5-6 


235 


E- 10-30 


53 


E-5-12 


399 


1-9-30 


73 


F-5-ie 


555 


E-8-30 


149 


F-5-74 


660 


E- 7—30 


£43 


F- 5-30 


684 


E-6-30 


399 


E-5-36 


654 


E-4-3U 


1167 


E-5-4.2 


537 


E- 3-30 


1842 


P-5-48 


392 


l -*-30 


2645 


r«s-54 


240 


E.-1-30 


3535 


F-5-60 


76.5 


A-5-50 


129 






a- o-5o 


316 






C-t-30 


469 






D-O-JC 


618 






F-w-iO 


699 






G-3-30 


617 






ts-O-30 


477 






1-0-30 


302 
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T abl e 10 • Hun 7 



Position 


Sfiturr ted 
act.ivlty(c/ffi) 


i osition 


8»turpt#d 

ectivlty(c/ra) 


L-13-30 


52 


J-o-50 


141 


F -10-30 


30 


E-5-0 


65 


-11-30 


61 


-5-6 


244 


F.- 10- 30 


78 


5 -5-12 


398 


E-9-30 


110 


-5-1 a 


531 


F-0-30 


193 


f -5- 4 


650 


-7-30 


268 


F-5-30 


702 


E-6- 50 


433 


E-5-3? 


689 


-4-50 


1094 


F-5-42 


541 


F -3-30 


1600 


F-5-48 


394 


r-I-30 


2330 


f- 5-54 


24 


F-1-.30 


5115 


'-5-50 


59 


A- 5-30 


136 






H-6-30 


315 






C-5-30 


492 






l -5-30 


562 






r-5-30 


723 






G-5-30 


644 






n-5-30 


4t;6 






1—5—30 


319 
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